Draft version February 25, 2013 

Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 



o 

(N 



Oh 

o 

in 



(N 

en 

o 

> 

• i— I 

X 



PROPERTIES OF SEQUENTIAL CHROMOSPHERIC BRIGHTENINGS AND ASSOCIATED FLARE RIBBONS 

Michael S. Kirk 1 ' 2,3 , K. S. Balasubramaniam 2-31 , Jason Jackiewicz 1 , R. T. James McAteer 1 , Ryan O. Milligan 4 

Draft version February 25, 2013 

ABSTRACT 

We report on the physical properties of solar sequential chromospheric brightcnings (SCBs) observed 
in conjunction with moderate-sized chromospheric flares with associated CMEs. To characterize these 
ephemeral events, we developed automated procedures to identify and track subsections (kernels) of 
solar flares and associated SCBs using high resolution Ha images. Following the algorithmic identifi- 
cation and a statistical analysis, we compare and find the following: SCBs are distinctly different from 
flare kernels in their temporal characteristics of intensity, Doppler structure, duration, and location 
properties. We demonstrate that flare ribbons are themselves made up of subsections exhibiting dif- 
fering characteristics. Flare kernels are measured to have a mean propagation speed of 0.2 km s and 
a maximum speed of 2.3 km s _1 over a mean distance of 5 x 10 3 km. Within the studied population 
of SCBs, different classes of characteristics are observed with coincident negative, positive, or both 
negative and positive Doppler shifts of a few km s~ x . The appearance of SCBs precede peak flare 
intensity by s=s 12 minutes and decay rj 1 hour later. They are also found to propagate laterally away 
from flare center in clusters at 41 km s _1 or 89 km s~ x . Given SCBs distinctive nature compared to 
flares, we suggest a different physical mechanism relating to their origin than the associated flare. We 
present a heuristic model of the origin of SCBs. 
Subject headings: Sun: chromosphere, Sun: coronal mass ejections (CMEs), Sun: flares 



1. INTRODUCTION 

The solar chromosphere exhibits three different classes 
of small scale intensity brightenings: flare-, plage-, and 
compact-brightenings. Although each is characterized 
by an enhanced temporal Ha brightness relative to a 
background quiet Sun, they each have distinct physical 
processes governing their spatial and temporal evolution. 
Typically brightenings have been identified and charac 
terized manually from a single data source (e.g. 



Kurt 



et al.||2000| |Ruzdjak et aL||1989| |Veronig et al7| 2002). 
However m order to form a better understanding of the 
underlying dynamics, data from multiple sources must 
be utilized and numerous similar features must be sta- 
tistically analyzed. 

This work focuses on flare brightenings and associated 
compact brightenings called sequential chromospheric 
brightenings (SCBs). SCBs were first observed in 2005 
and appear as a series of spatially separated points tha t 
brighten in sequence ( Balasubramaniam et al. 2005). 
SCBs are observed as multiple trains of brightenings m 
association with a large-scale eruption in the chromo- 
sphere or corona and are interpreted as progressive prop- 
agating disturbances. The loci of brightenings emerge 
predominantly along the axis of the flare ribbons. SCBs 
are correlated with the dynamics which cause solar flares, 
coronal restructuring of magnetic fields, halo C MEs, EIT 
waves, and chromosphe ri c sympathetic fl a ring ^Balasub- 
ramaniam et al.||2005|). iPevtsov et al.l (120076 demon- 
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strate that SCBs have properties consistent with aspects 
of chromospheric evaporation. 

This article presents a new description of the dynam- 
ical properties of S CBs resulting fro m applying a new 
automated method ( |Kirk et al.| 2011 ) of identifying and 
tracking SCBs and associated flare ribbons. This track- 
ing technique differs from previous flare tracking algo- 
rithms in that it identifies and tracks spatial and tempo- 
ral subsections of the flare and all related brightcnings 
from pre-flare stage, through the impulsive brightening 
stage, and into their decay. Such an automated mea- 
surement allows for tracking dynamical changes in in- 
tensity, position, and derived Doppler velocities of each 
individual subsection. The tracking algorithm is also 
adapted to follow the temporal evolution of ephemeral 
SCBs that appear with the flare. In Section|2]we describe 
the data used to train the algorithm and the image pro- 
cessing involved in the detection routine. In Section [3] we 
present the results of tracking the evolution of flare Ker- 
nels through an erupting flare. In Section [4] we present 
the application of the tracking algorithm to ephemeral 
SCBs. We present a physical interpretation of the SCBs 
and provide a heuristic a model of the origin of SCBs in 
Section[5] Finally, in Section[6]we discuss implications of 
these results and provide future direction for this work. 

2. DATA AND DATA PROCESSING 

In this study we use chromospheric Ha (6562.8 A) im- 
ages from the US AF's Improved Sol ar Observing Optical 
Network (ISOON; lNeidig et al.|1998[ ) prototype telescope 
to study flare ribbons and SCBs. ISOON is an auto- 
mated telescope producing 2048x2048 pixel full-disk im- 
ages at a one-minute cadence. Each image has a 1.1 arc- 
second spatial sampling, is normalized to the quiet Sun, 
and corrected for atmospheric refraction (Figures fTfa), 
[2]ja), and|3^a)). Nearly coincident to the spectralline 



Table 1 

The events selected for study. Each event was visually identified 

as having a two-ribbon configuration and SCBs associated with 

the flare. 



Date Start (UT) Duration (h) 
9 Nov. 2004 16:59 05 

6 May 2005 16:03 2.1 

13 May 2005 16:13 1.3 



Flare Class CME 

M8.9 Halo 

C8.5 Halo 

M8.0 Halo 



center images, ISOON also records Ha ±0.4 A off-band 
Doppler images, from which Dopplcr signals are derived. 
For this study, we chose to apply the brightening detec- 
tion a l gorith ms to three flares where Balasubramaniam 



et aL (20061 had previously identified SCB events (Ta- 
ble HI . Each of the events selected has a two-ribbon con- 
figuration and an associated halo CME. Of these, both 
the 6 May 2005 and 13 May 2005 events were both lo- 
cated near disk center while the 9 November 2004 event 
was near the western limb. Images were extracted from 
the archive from ±3.5 hours of the eruption start time, 
yielding a data cube with sa 400 images for each event. 
An Ha relative intensity curve, an Ha maximum inten- 
sity curve, and Geostationary Operational Environmental 
Satellites (GOES) hard and soft x-ray fluxes are plotted 
to characterize the flare (Figures [IJd) , [2JJd), and[3]jd)) 
as described in Kirk et al. ( 2011[ ). 

Each ISOON image is reconditioned to remove so- 
lar limb-darkening. The images are then de-projected 
into conformal coordinates using a Guyou project i on (an 
oblique aspect of the Peirce projection, Peirce 1879), 
which removes the projection effects of imaging the so- 
lar sphere. Each image is then cropped to the region of 
interest (ROI) (Figures jljb), gb), and^b)). The set 
of images are aligned using a cross-correlation algorithm 
eliminating the rotation effects of the Sun. Frames that 
contain bad pixels or excessive cloud interference are re- 
moved. 

The red and blue wings of ISOON Doppler images are 
each preprocessed using the same technique as applied 
to the line center images described above (Figures |lj c ) , 
2Tc), andpjlc)). In order to produce velocity measure- 
ments, a Doppler cancelation technique is employed in 
which red images are su btracted from blue images (for 
a modern example see: Connes 19851. Typically, red 
and blue images are separated by about 4 seconds and 
no more than 5 seconds. A mean zero redshift in the 
subtracted Doppler image in the quiet-Sun serves as a 
reference. In dynamical situations such as solar flares, 
Ha profiles are often asymmetric especially in emission 
making this Doppler technique invalid. To avoid asym- 
metric profiles, only values derived outside of the flaring 
region are considered such that the areas of interest at 
±0.4 A are still in absorption e ven if raised in intensity 
(Balasubramaniam et al. 2004). The values from the 
subtracted images are translated into units of km s _1 by 
comparison to a measured response in spectral line shifts, 
which are calibrated against an intensity difference for 
the ISOON telescope. In this context, we interpret the 
Doppler shift as a line of sight velocity. Thus, Doppler 
velocity [vrj] is defined as: 



2( 



VD — k (Ibluo — Zred) . 



(1) 



where I is the measured intensity and k is a linear fitted 
factor that assumes the intensity changes due to shifts in 



the symmetric spectral line, which can be attributed to 
a Doppler measure as a first-order approximation. The 
linear factor [k] has been independently determined for 
the ISOON telescope by measuring the full H a spectral 



line Doppler shift across the entire solar disk (Balasub- 
|ramaniam|201lj ). 

2.1. Detection and tracking 

An animated time series of sequential images cover- 
ing an erupting flare reveals several physical characteris- 
tics of evolving ribbons: the ribbons separate, brighten, 
and change their morphology. Adjacent to the eruption, 
SCBs can be observed b rightening a n d d imming in the 



vicinity of the ribbons. Kirk et al. (20111 describe in 



detail techniques and methods used to extract quanti 
ties of interest such as location, velocity, and intensity 
of flare ribbons and SCBs. The thresholding, image en- 
hancement, and feature identification are tuned to the 
ISOON data. The detection and tracking algorithms are 
specialized for each feature of interest and requires phys- 
ical knowledge (e.g. size, peak intensity, and longevity) 
of that feature being detected to isolate the substructure. 

Briefly, the detection algorithm first identifies candi- 
date bright kernels in a set of images. In this context, we 
define a kernel to be a locus of pixels that are associated 
with each other through increased intensity as compared 
with the immediately surrounding pixels. Each kernel 
has a local maximum, must be separated from another 
kernel by at least one pixel, and does not have any prede- 
termined size or shape. Next, the algorithm links time- 
resolved kernels between frames into trajectories. A filter 
is applied to eliminate inconsistent or otherwise peculiar 
detections. Finally, characteristics of bright kernels are 
extracted by overlaying the trajectories over complemen- 
tary datasets. To aid in this detection, tracking software 
developed by Crocker, Grier and Weeks was used as a 
foundation and modifie d to fit the needs of this projeclQ 
( |Crocker fc Gner1p96l ). 

In order to characterize a kernel, we calculate its inte- 
grated intensity, radius of gyration, and eccentricity. The 
eccentricity of the kernel (as defined by its semi-major [a] 
and semi-minor [b] axis) is calculated using, 
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The integrated intensity of a given kernel is defined by 

where A(x + i,y + j) is the intensity of the pixel located 
at (x + i, y + j), the kernel has a brightness weighted 
centroid w ith coordinates (x.y), a nd v is the radius of 
the mask (Crocker & Grier 1996). The radius of the 
mask is chosen to be eight pixels m this case. A radius 
of gyration, R g , is related to the moment of inertia [7] 



by using: 



I = ^2m k rl = mR 2 g , 

k 

r = (a + b)/2, 



(4) 



(5) 



1 Crocker's software is available online at www. physics. emory. 
edu/$\sim$weeks/idl/ 



where mk is the mass of particle k, m is the total mass 
of the system, and r is the distance to the rotation axis. 
In this context, we interpret the radius of gyration to be 



R 2 g = 
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where m is the integrated int ensity as defined in Equa- 
tion [3] ( Crocker & Gricr 1996). Figure [7] diagrams a hy- 
pothetical kernel. The results of this algorithm identify 
and characterize several kernels in the flaring region (Fig- 
ure Eft. The total number of flare kernels detected is typi- 
callybetween 100 - 200, while the number of SCB kernels 
detected is typically three to four times that number. 

3. FLARE RIBBON PROPERTIES 

The properties of kernels identified in flare ribbons can 
be examined in two ways: (i) each kernel can be consid- 
ered as an independent aggregation of compact bright- 
enings or (ii) kernels can be considered as dependent on 
each other as fragments of a dynamic system. Each cate- 
gory brings about contrasting properties. If one considers 
the kernels as independent elements, this provides a way 
to examine chan ges t o subsections of the flare and is dis- 
cussed in Section [XT] Associating kernels with their con- 
textual surroundings allows a way to examine the total 
evolution of the flare without concern of how individual 
kernels behave. This type of examination is addressed in 
Section l3~2l 

3.1. Qualities of individual kernels 

A sample of six kernels from the 13 May 2005 event, 
letters U - X in Figure [5] is representative of the ma- 
jority of flare kernels tracked. The kernels were selected 
from three different regions of this flare. The normalized 
Ha intensity of each of these kernels is shown in Fig- 
ure [6j In contrast with the integrated flare light curve, 
an individual flare kernel often has a shorter lifetime. 
This is because an individual kernel may first appear in 
the impulsive phase of the flare (as demonstrated in ker- 
nels W - Y) while other kernels may disappear as the 
flare begins to decay (as in kernels U and Z). Most likely 
these kernels merged with one another, at which point 
their unique identity was lost. Kernels W and X are the 
only two that show the characteristic exponential dim- 
ming found in the reference curves in Figuren](d). All of 
these kernels have peak intensities within a few minutes 
of the peak of the total flare intensity and have sustained 
brightening above background levels for up to an hour. 

The integrated speed of displaced kernels provides con- 
text to the evolution of intensities. We define the inte- 
grated speed of displacement to be the sum over all time 
steps of the measured velocity of an individual kernel and 
thus small velocity perturbations are minimized by the 
sum. The peak integrated speed measured for each ker- 
nel peaks at ~ 2.3 km s _1 and has a mean of ~ 0.2 km 
s _1 . There is significant motion along the flare ribbons 
as well as outflow away from flare center. The motions 
are complex but generally diverge. The total spatial (an- 
gular) displacement of each flare kernel peaks at ~ 20.7 
Mm and has a mean of ~ 5.0 Mm. Flare kernels that ex- 
hibited the greatest speeds did not necessarily have the 
greatest displacements. As the flare develops, most of the 
motion in flare ribbons indicated is synchronous with the 



separation of the ribbons. A few tracked kernels indicate 
motion along flare ribbons. These flows along flare rib- 
bons are consistently observed in all three flares. 

3.2. Derived flare quantities 

Examining trajectories of flare kernels offers insight 
into the motions of the flare ribbons as they evolve 
through the eruption. Figure [HI shows a time series of 
images centered about the peak-time of the flare. The 
velocities of the detected flare kernels are superposed as 
vectors. Notice the initial outflow of the two ribbons near 
the flare peak at 16:49 UT. But, as the flare continues to 
evolve, there is more motion along the flare ribbons. This 
is probably a result of the over-arching loops readjusting 
after the reconnection event. Beneath each image is a 
histogram of the distribution of velocities at each time 
step. The mean apparent lateral velocity of each ker- 
nel remains below 0.5 km s _1 throughout the flare even 
though the total velocity changes significantly. Qualita- 
tively, the bulk of the apparent motion appears after the 
flare peak intensity at 16:49 UT. 

Grouping the kernels' velocities into appropriate bins 
provides another way to analyze the dynamics of the flare 
(Figure |9|. The kernels with apparent lateral speeds 
above 0.4 km s _1 are highly coincident with the peak 
intensity of the flare. The number of these fast moving 
kernels peaks within a couple minutes of flare-peak, and 
then quickly decays back to quiet levels. The velocity 
bins between 0.1 and 0.4 km s _1 show a different sub- 
structure. These velocity bins peak rs 30 minutes after 
the peak of the flare intensity. They show a much slower 
decay rate, staying above the pre-flare velocity measured 
throughout the decay phase of the Ha flare curve. 

Figure 10 (left column) shows the integrated unsigned 
kernel velocity (speed) for each image. In general, the 
evolution of the integrated kernel speed has a shape sim- 
ilar to the intensity curve. For the 6 May 2005 event, the 
speed curve is more similar to the GOES x-ray intensity 
curve than the Ha curve. For the November 2004 event, 
the limb geometry of the flare, as well as the noise in the 
original dataset, results in a noisy curve. Despite these 
difficulties, a clear increase in total speed is apparent near 
the peak of the flare. A continuum integrated speed level 
of 1-2 km s _1 is consistent in each event, indicating the 
quality of data and of the tracking software. The peak 
integrated speed is about 15 km s" 1 for all three flares. 

Examining the change in the speed of each kernel be- 
tween successive images yields a derived acceleration for 
each tracked kernel. The integrated unsigned accelera- 
tion is plotted in Figure 10 (middle column). Both the 13 
May 2005 and 6 May 2005 events show peak acceleration 
after the peak of the intensity curve. Hence the major- 
ity of the acceleration in the apparent motion of the flare 
comes in the decay phase of Ha emission, in concert with 
the formation of post-peak flare loops, as seen in movies 
of TRACE images of other flares. The peak acceleration 
appears uncorrelated with the strength of the flare since 
the peak of the M class 13 May event is just over 300 km 
s~ 2 , while the peak of the C class flare on 6 May is nearly 
400 km s~ 2 . The acceleration curve for the 9 November 
flare has too much noise associated with it to decisively 
determine the peak value. 



The right column of Figure 10 shows a derived kinetic 
energy associated with the measured motion of the flare 



kernels. This was accomplished by assuming a chromo- 
spheric density [a] of 10 -5 kg m -3 and a depth [h] of 
1000 km. The mass contained in each volume is then 
calculated by 

m = TrhaR 2 g y/l-e 2 (7) 

where R g is defined in Equation [6] and e is defined in 
Equation [2j The kinetic energy under a kernel is defined 
in the standard way to be 



El 



(8) 



where v (speed) is a measured quantity for each kernel. 
The derived kinetic energy curves show a ~ 30 fold in- 
crease in kinetic energy during the flare with a decay 
time similar to the decay rate of the Ha intensity. The 
measure of a kernel's kinetic energy is imprecise because 
the measured motion is apparent motion of the underly- 
ing plas ma. In the model of two ribbon solar flares put 



forth by Priest & Forbes ([2002), the apparent velocity 
is a better indicator ot rates or magnetic reconnection 
rather than of plasma motion. Despite this caveat, the 
derived kernel kinetic energy is a useful measure because 
it combines the size of the kernels with their apparent 
velocities to characterize the flare's evolution, thus rep- 
resenting the two in a quantitative way. 

The non-thermal 25-50 keV emission between 16:42- 
16:43 UT from the Reuven Ra maty High Energ y Solar 



Spectroscopic Imager (RHESSI, Lin et al. 2002 ) of the 



13 May flare is contoured in Figure [TT] over an Ha image 
taken at 16:42 UT. The high-energy emission is centered 
over the flare ribbons but is discontinuous across the rib- 
bon. There are three localized points from which the ma- 
jority of the x-ray emission comes. The parts of the flare 
ribbons exhibiting the most displacement show generally 
lower x-ray intensity than the more stationary segments 
of the flare ribbons. 

4. SCB PROPERTIES 

The measured properties of SCBs can be studied using 
two scenarios: each SCB kernel is considered as indepen- 
dent and isolated compact brightcnings; each compact 
brightening is considered to be fragments of a dynamic 
system that has larger structure. This is similar to the 
approach used for flares. Again, considering both sce- 
narios has its benefits. Identifying kernels as indepen- 
dent elements reveals structural differe nces in compact 
brightenings and is discussed in Section |4~Tj Associating 
kernels with their contextual surroundings provides some 
insight into the physical structure causing SCBs. An ag- 



gregate assessment of SCBs is addressed in Section 4.2 



To mitigate the effects of false or marginal detections, 
only SCBs with intensities two standard deviations above 
the mean background intensity are considered for char- 
acterizing the nature of the ephemeral brightening. 

4.1. Qualities of individual SCBs 

SCBs, although related to erupting flare ribbons, are 
distinctly different from the flare kernels discussed in Sec- 
tion [3j Six SCBs are chosen from the 13 May 2005 event 
as an example of these ephemeral phenomena. The lo- 
cations of these six events are letters A - F in Figure [5| 
The Ha normalized intensity of each of these kernels is 



shown in Figure [6] The SCB intensity curve is signif- 
icantly different from the flare kernel curves shown on 
the left side of the figure. SCB curves are impulsive; 
they have a sharp peak and then return to background 
intensity in the span of about 12 minutes. Nearly all of 
the SCBs shown here peak in intensity before the peak of 
the flare intensity curve, shown as a vertical dashed line 
in Figure [6} SCBs B and E both appear to have more in- 
ternal structure than the other SCBs and last noticeably 
longer. This is most likely caused by several unresolved 
SCB events occurring in succession. 

Nearly all SCBs have Ha line center intensities 40 - 
60% ab ove the mean intensity of the quiet chromosphere 
(Figure 12). From this increased initial intensity, SCBs 
brighten to intensities 75 - 130% above the quiet Sun. 

Examining Doppler velocity measurements from the 
SCB locations reveals three distinct types of SCBs (Fig- 



12 ). A type la SCB has an impulsive intensity profile 



and an impulsive negative Doppler profile that occurs si- 
multaneously or a few minutes after the peak brighten- 
ing. In this study, a negative velocity is associated with 
motion away from the observer and into the Sun. A type 
lb SCB has a similar intensity and Doppler profile as a 
type la but the timing of the impulsive negative velocity 
occurs several minutes before the pea k in tensity of the 
SCB. The type lb SCB shown in Figure 12 has a negative 
velocity that peaks 10 minutes before the peak intensity 
and returns to a stationary state before the Ha line cen- 
ter intensity has decayed. A type II SCB has a positive 
Doppler shift perturbation that often lasts longer than 
the emission in the Ha intensity profile. The timing of 
both are nearly coincidental. A type III SCB demon- 
strates variable dynamics. It has a broad Ha intensity 
line center with significant substructure. A type III SCB 
begins with a negative Doppler profile much like a type 
I. Before the negative velocity perturbation can decay 
back to continuum levels, there is a dramatic positive 
velocity shift within three minutes with an associated 
line center brightening. In all types of SCBs the typical 
magnitude of Doppler velocity perturbation is between 
2 and 5 km s _1 in either direction perpendicular to the 
solar surface. 

Of all off-flare compact brightenings detected using the 
automated techniques, 59% of SCBs have no discernible 
Doppler velocity in the 6 May event and 21% in the 13 
May event. Out of the SCBs that do have an associated 
Doppler velocity, the 6 May event has 35% of type I, 54% 
of type II, and 11% of type III. The 13 May event has 
52% of type I, 41% of type II, and 17% of type III. When 
totaling all of the Doppler velocity associated SCBs be- 
tween both events, 41% are of type I, 46% of type II, 
and 13% of type III. Both line center brightenings and 
Doppler velocities are filtered such that positive detec- 
tions have at least a two standard deviation peak above 
the b ackground noise determined in the candidate detec- 
tion ( |Kirket al.|2011 ). 



4.2. SCBs in aggregate 

Examining SCBs as a total population, the intensity 
brightenings have a median duration of 3.1 minutes and 
a mean duration is 5.7 minutes (Figure 13). The dura- 
tion is characterized by the full-width halt-maximum of 
the SCB intensity curve (examples of these curves are 
shown in Figure H3J) . A histogram of the distribution of 



the number of SCB events as a function of duration shows 
an exponential decline in the number of SCB events be- 
tween 2 and 30 minutes (Figure 13 1. The duration of 



SCBs is uncorrelated with both distance from flare cen- 
ter and the peak intensity of the SCB. 

Examining individual SCBs' distance from the flare 
center as a function of time provides a way to extract 
the propagation trends of SCBs around the flare (Fig- 



ure 



14|. The bulk SCBs appear between 1.2 and 2.5 Mm 
away from flare center in a 1 hour time window. At the 
extremes, SCBs are observed at distances up to « 5 Mm 
and several hours after the flare intensity maximum. On 
the top two panels of Figure [Ml the shade of the mark 
corresponds to the intensity ofthe SCB measured where 
the lighter the mark, the brighter the SCB. The center 
of the flare is determined by first retaining pixels in the 
normalized dataset above 1.35. Then all of the images 
are co-added and the center of mass of that co-added 
image is set to be the "flare center." The SCBs tend to 
clump together in the time-distance plot in both the 6 
May and 13 May events. Generally, the brighter SCBs 
are physically closer to the flare and temporally occur 
closer to the flare peak. This intensity correlation is weak 
and qualitatively related to distance rather than time of 
brightening. Statistics from the November 9, 2004 event 
are dominated by noise. 

Identifying SCBs with negative Doppler shifts within 
3 minutes of the Ha peak intensity limits SCBs to type 
I only. With this method of filtering, some trends be- 
co mc m ore apparent in a time-distance plot of SCBs (Fig- 
ure [14 bottom panels). In this view of SCBs, they tend 



to cluster together in time as well as distance, with two 
larger groups dominating the plots. 

To fit a slope to these populations, a forward-fitting 
technique is employed similar to a linear discriminate 
analysis. This method requires the user to identify the 
number of groups to be fitted and the location of each 
group. The fitting routine then searches all linear combi- 
nations of features for the next "best point" to minimize 
the chi-square to a regression fit of the candidate group. 
Repeating this method over all the points in the set pro- 
duces an ordered set of points that when fit, have an 
increasing chi-square value. A threshold is taken where 
the derivative of the chi-square curve increases to beyond 
one standard deviation, which has the effect of identify- 
ing where the chi-square begins to increase dramatically. 
Running this routine several times minimizes the effect 
of the user and provides an estaminet of the error asso- 
ciated with the fitted line. This method has two caveats. 
First, this fitting method relies on the detections having 
Poissonian noise. This is not necessarily correct since 
the detection process of compact brightenings introduces 
a selection bias. Second, the fitting method makes the 
assumption that no acceleration occurs in the propaga- 
tion of SCBs. This is a reasonable approximation but 
from studies of Moreton wave s in the chromosphere (e.g. 
Balasubramaniam et al.|2010 ), a constant velocity prop- 
agation is unlikely. 

Applying the forward-fitting technique these two data 
sets yields two propagation speeds: a fast and a slow 
group. The 6 May event has propagation speeds of 71 ±2 
km s^ 1 and 49 ± 2 km s _1 . The 13 May event has 
propagation speeds of 106 ± 7 km s~ x and 33 ± 4 km 
s _1 . Beyond these detections, the 13 May event has a 



third ambiguous detection with a propagation speed of 
— 7 ± 5 km s _1 . The implications of these statistics are 
discussed in Section [6] 

The x-ray intensity measured with RHESSI tempo- 
rally covers a part of the flaring event on 13 May 2005. 
RHESSI data coverage begins at 16:37 UT during the 
impulsive phase of the flare. Comparing the integrated 
x-ray intensity curve between 25 - 50 keV to the ag- 
gregate of SCB intensities integrated over each minute 
yields some similarities (Figure 15). The x-ray intensity 
peaks about a minute after the integrated SCBs reach 
their maximum intensity. The decay of both the x-ray 
intensity and SCB integrated curve occur approximately 
on the same time-scale of « 50 minutes. 

5. PHYSICAL INTERPRETATION OF RESULTS 

From the results of the tracking method, we can gain 
insight into the physical processes associated with SCBs 
and flares. First, the flare ribbons studied are both spa- 
tially and temporally clumped into discrete kernels. The 
relatively smooth ribbon motions observed in Ha and 
the nearly ideal exponential decay characterized in x-ray 
flux measurements belie the substructure that mak es up 
a two- ribbon flare. This is discussed in Section IBTTl Sec- 
ond, the brightening in the flare ribbons is caused by a 
distinctly different physical process than the one causing 
SCBs. In this sense, SCBs are different than micro-flares. 
Our interpretatio n of the physical origin of SCBs is ad- 
dressed in Section [521 



5.1. Two ribbon chromo spheric flares 

Flare kernels are observed to appear and disappear as 
the unde rlying flare ribbons evolve, as discussed in Sec- 
tion [XT] Examining individual kernel structure suggests 
there is substructure within a flare ribbon whose ele- 
ments impulsively brighten and dim within the bright- 
ness that encompasses the intensity curve. These results 
support the premise that flares are m ade up of several 
magn etic field lines reconnecting (e.g., |Priest fc Forbes| 
2002 1 . There is no evidence to claim an individual Hare 
kernel is directly tracking a loop footpoint. Within a 
tracked flare kernel, multiple coronal reconnection events 
are superimposed to produce the observed asymmetries 
a single flare kernel's light cur ve. 



Maurya & Amba stha| (2010| tracked subsections ofthe 
28 October 2003 XI 7 flare as it evolved. They reported 
peak speeds ranging from ~ 10 — 60 km s _1 (depending 
on the spatially tracked p art of the flare) over the ob- 
served span of 13 minutes. |Maurya fc Ambastha (20101 
also reported that the total apparent distances the rib- 
bons traveled were ~ 10 4 km. The 28 October flare is 
several orders of magnitude greater in GOES x-ray in- 
tensity than the flares considered for this study. In the 
present study, despite this difference, peak speeds of flare 
kernels observed in this work are measured at ~ 2.3 km 
s" 1 and the mean velocity of all flare kernels is ~ 0.2 
km s _1 . The maximum distance flare kernels traversed 
was ~2x 10 4 km and the average distance traveled was 
~5x 10 3 km; both of these values are similar to the 28 
October 2003 event. Although the velocities of the flare 
ribbons studied here are at least an order of magnitude 
less, the two to three orders of magnitude difference in 
the GOES peak intensity between the 28 October flare 
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and the flares studied here implies that the apparent ve- 
locities and distance traveled by flare ribbons do not scale 
linearly with the strength of the eruption. 

The bulk of the apparent motion and acceleration in 
the flare ribbons are observed after the peak of the flare 
intensity. Integra ting this into the dynamic al model of a 
two- ribbon flare (Demoulin & Priest 1988) implies that 
the peak intensity of the flare occurs m the low lying 
arcade and loses intensity as the x- point reconnection 
progresses vertically to higher levels ([Baker et al.||2009 |. 
The measured divergence of the ribbons dominate the 
motion but there is significant apparent flow tangential 
to the flare ribbons. This predicates that there is lateral 
propagation to the x-point as well as vertical propaga- 
tion. 



5.2. A physical model of SCBs 

Figure [16] is a proposed phenomenological model of the 
overlying physical topology. SCBs are hypothesized to be 
caused by electron beam heating confined by magnetic 
loop lines over-arching flare ribbons. These over-arching 
loops are analogous to the higher lyin g, unsheared tether s 
in the breakout model of CMEs ( |Antiochos et al.| [l999). 
As the flare erupts, magnetic reconnection begins m a 
coronal x-point. The CME escapes into interplanetary 
space, the remaining loop arcade produces a two ribbon 
flare, and the tethers reconnect to a new equilibrium 
position. The tether reconnection accelerates trapped 
plasma which impacts the denser chromosphere causing 
observed brightening. This description implies the driver 
of SCBs is t he eruption of an CME. Balasubramaniam 
et al. (2006) found a strong correlation between CMEs 
and the presence of SCBs. 

A simple loop configuration without localized diffusion 
or anomalous resistance implies that the length of the 
loop directly proportional to the travel time of the elec- 
tron beam. This means t hat the different propagation 
groups observed in Figure 14 result from different phys- 
ical orientations of over-arching loops. The time it takes 
to observe a chromospheric brightening can therefore be 
described as: 

L 
tscb=7 — i (9) 



where tscb is the time it takes for accelerated plasma to 
travel from the coronal x-point, along the magnetic loop 
lines, impact the chromosphere, and produce a bright- 
ening, L is the length of the tether before reconnection, 
v e is the electron velocity along the loop line, and 7 is 
a function of changing diffusion and resistance between 
chromosphere and corona. In physical scenarios, v e is 
likely a combination of the Alfven speed and the elec- 
tron thermal velocity. 

The slopes plotted in Figure [14] show two different 
propagation speeds for SCBs: 33—49 km s~ x and 71 — 106 
km s . The sound speed in the chromosphere is ap 



proximately c s ~ 10 km s (Nagashima et al. 2009), 
while the Alfven speed in the upper chromosphere is ap- 
proximated to b e between va ~ 10 — 100 km s^ 1 (As 



chwandcn 2005). Both of the SCB propagation speeds 
are significantly above the approximated sound speed, 
however they fall reasonably well into the range of the 
Alfven speed. Assuming the propagation speed of non- 
thermal plasma along the loop line is consistent between 



flare loops, the different propagation speeds would then 
imply different populations of loops being heated as the 
flare erupts. 

Examining single SCBs, the coincident Doppler re- 
coil with the Ha intensity presents a contradiction. If 
an SCB is a n example of compact chro mospheric evap- 
oration (e.g. Dennis & Schwartz 1989), then the only 
Doppler motion should be outward, the opposite of what 
is observed. Figure[l7]presents a possible solution to this. 
Since the scatteringlength of electrons is significantly 
smaller than that of protons, the electron beam impacts 
the mid chromosphere and deposit energy into the sur- 
rounding plasma while the protons penetrate deeper. 
This deposited heat cannot dissipate effectively through 
conduction or radiation and thus expands upward into 
the flux tube. To achieve an expansion with rates of a 
few km s _1 as observed, the chromospheric heating rate 



must be below Eh < 10 erg cm s (Fisher et al. 



1984). As a reaction to this expansion, a reaction wave 



propagates toward the solar surface in the opposite direc- 
tion of the ejected plasma. Since the Doppler measure- 
ments are made in the wings of the Ha line, the location 
of the Doppler measurement is physically closer to the 
photosphere than the Ha line center. Thus the observer 
sees the heating of the Ha line center and coincidentally 
observes the recoil in the lower chromosphere. 

In type II SCBs an up-flow is observed. This is an 
example of a classic chromospheric evaporation where 
the heated plasma in the bulk of the chromosphere is 
heated and ablated back up the flux tube. In contrast, 
the type III class of SCB shown in Figure [12] present 
an interesting anomaly to both the other observed SCBs 
and the model proposed in Figure 17 Since there is an 
initial down-flow, the beginning state of type III SCBs 
are similar to type Is. As the recoil is propagated, a 
continual bombardment of excited plasma (both protons 
and electrons) impacts the lower chromosphere, causing 
ablation and changing the direction of flow. 

6. SUMMARY AND CONCLUSIONS 

Two conclusions can be made about the flares stud- 
ied in this work: first, the asymmetrical motions of the 
flare ribbons imply the peak flare energy occurs in the 
low lying arcade, and second, flare related SCBs appear 
at distances on the order of 10 5 km and have pro perties 
of sites of compact chro mospheric evaporation 
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rya & Ambastha (20101 estimated reconnection rates 
from a measurement ot the photospheric magnetic field 
and the apparent velocity of chromospheric flare ribbons. 
Therefore, it would be possible to estimate the recon- 
nection rates using this technique with the addition of 
photospheric magnetograms. Associating vector mag- 
netograms with this technique and a careful consider- 
ation of the observed Doppler motions underneath the 
ribbons would provide a full 3D method for estimating 
the Lorentz force for subsections of a flare ribbon. 

SCBs originate during the impulsive rise phase of the 
flare and often precede the Ha flare peak. SCBs are 
found to appear with similar qualities as comp act chro- 
mospheric a blation confirming the results of IPevtsov 



et al. (2007). The heuristic model proposed in Figure 17 
requires unipolar magnetic field underneath SCBs. The 
integration of high resolution magnetograms and a field 
extrapolation into the chromosphere would confirm this 



chromospheric evaporation model. One consequence of 
chromospheric evaporation is that the brightening should 
also be visible in EUV and x-ray observations due to non- 
thermal protons interacting with the lower chromosphere 
or photosphere. 

SCBs are a special case of chromospheric compact 
brightening that occur in conjunction with flares. The 
distinct nature of SCBs arises from their impulsive 
brightenings, unique Doppler velocity profiles, and origin 
in the impulsive phase of flare eruption. These facts com- 
bined demonstrate that SCBs have a non-localized area 
of influence and are indicative of the conditions of the en- 
tire flaring region. They can possibly be understood by 
a mechanism in which a destabilized overlying magnetic 
arcade accelerates electrons along magnetic tubes that 
impact a denser chromosphere to result in an SCB. This 
distinguishes SCBs from the flare with which they are 
associated. In a future work, we plan to address physical 
mechanisms to explain the energetic differences between 
SCBs and flares to explain the coupled phenomena. 

The authors thank: (1) USAF/AFRL Grant FA9453- 
11-1-0259, (2) NSO/AURA for the use of their Sunspot, 
NM facilities, (3) AFRL/RVBXS, and (4) Crocker and 
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Figure 1. The two-ribbon flaring event from 13 May 2005. A shows an example of a calibrated Ha ISOON image with the region of 
interest (ROI) highlighted. The ISOON images are 2048 X 2048 pixels with each pixel having a 1.1 arcsec resolution. B is the ROI after 
preprocessing. This ROI coves 576 X 576 pixels corresponding to 2.12 X 10 11 km 2 on the solar surface. C is a Doppler measure of the same 
ROI as part B. The Doppler velocity image is created out of ISOON Ha off-band images described in Section pjwhere the velocity ranges 
from -26.6 to 21.5 km s from black to white. D shows intensity curves over the time period of interest: 13:00-20:02 UT as described in 
detail in the text. 
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Figure 2. The two-ribbon flaring event from 6 May 2005. A shows an example of a calibrated Ha ISOON image with the region of 
interest (ROI) highlighted. The ISOON images are 2048 X 2048 pixels with each pixel having a 1.1 arcsec resolution. B is the ROI after 
preprocessing. This ROI coves 768 X 768 pixels corresponding to 3.77 X 10 11 km 2 on the solar surface. C is a Doppler measure of the same 
ROI as part B. The Doppler velocity image is created out of ISOON Ha off-band images described in Section [2] where the velocity ranges 
from -12.5 to 6.9 km s — 1 from black to white. D shows intensity curves over the time period of interest: 12:28^22:47 UT as described in 
detail in the text. 
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Figure 3. The two-ribbon flaring event from 9 November, 2004. A shows an example of a calibrated Ha ISOON image with the region of 
interest (ROI) highlighted. The ISOON images are 2048 X 2048 pixels with each pixel having a 1.1 arcsec resolution. B is the ROI after 
preprocessing. This ROI coves 600 X 600 pixels corresponding to 2.30 X 10 11 km 2 on the solar surface. C is a Doppler measure of the same 
ROI as part B. The Doppler velocity image is created out of ISOON Ha off-band images described in Section [2] where the velocity ranges 
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Figure 4. Example ISOON Ha images from the 13 May 2005 two ribbon flare. The detection algorithm is applied to the image series to 
extract either ribbon brightening or compact brightening. The left image shows one time slice of individual SCBs identified external to the 
flare ribbon. The right image shows the same time slice with brightcnings identified along the flare ribbon. Each detection is made up of 
a small matrix of pixels associated with the brightening called a kernel. 
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Figure 5. An image from the 13 May 2005 event. Letters A - F point to the locations several SCB kernels that are shown in more detail 
in Figure [6] Letters U - Z point to the locations of a sample of several flare kernels that are also shown in more detail in Figure [6] Note 
that since nare kernels move with the evolving flare ribbon, the locations of these kernels are only accurate for this image. In contrast, 
SCB kernels are stationary and remain in the same location in any image. 
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Figure 6. Ha normalized intensity curves for each of the kernels indicated in Figurepjlby the associated letters. On the left, the tracked 
flare kernels are shown. The normalized intensity curves for the SCB kernels, shown on the right, are stationary and therefore tracked 
throughout the flaring window. In both panels, the vertical dashed line indicates the time at which the flare intensity peaks (16:49 UT) 
derived from integrating all flare kernel intensities at each time step. 
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Figure 7. A diagram of a hypothetical kernel (either flare or SCB) used for classifications in this paper. Each pixel is labeled with a 
representative intensity [Ai\. The semi-major axis [a] and semi-minor axis [b] are shown. This kernel has the following properties: total 
intensity, mo = 43, semi-major axis, a = 4, semi-minor axis, b = 2, eccentricity, e = 0.87, centroid, (x,y) = (0.26,0.18), and radius of 
gyration, R g = 1.9. The maximum intensity, in this case 5, is defined as having pixel coordinates (0,0). 
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Figure 8. A time series of images with the apparent lateral velocity of each tracked flare kernel at that time step plotted over the Ho 
flare ribbons. A histogram of the number of flare kernels as a function of velocity corresponding to each time step is plotted underneath. 
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Figure 9. The number of elements in each apparent lateral velocity bin plotted as a function of time. The vertical dashed line marks the 
time of the peak flare intensity. The velocity bins are 0.1 km s _1 in width. The smallest bin, 0-0.1 km s — 1 , is excluded because of a base 
level of noise inherent in the detection routine. 
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Figure 10. The integrated unsigned propagation speed, acceleration, and energy of flare kernels as a function of time for each of the 
three events studied. The vertical dashed line indicates the maximum flare intensity for the flaring event. The first column shows the 
integrated speed curves. The second column shows the acceleration of the flare ribbons. The third column shows the integrated kinetic 
energy, equation [8] as a function of time for each of the three events studied. The 2004 flare had significant noise in the original data and 
a location on the umb contributing to the inconsistent measurements. 
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Figure 11. An Ha image of the May 13, 2005 flare at peak intensity is shown with the reconstructed RHESSI 25 — 50 kcV pass band 
integrated between 16:42-16:43 UT contoured on top at 10, 30, 50, 70, and 90% of the maximum intensity. The field of view of the RHESSI 
image [-220:-95, 180:305] is confined within the flaring region. 
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Figure 12. The line center intensity and Doppler velocity measurements for three different types of SCBs observed in both 2005 flares. 
The vertical dashed line marks the peak Ho intensity of the associated SCB. A negative Doppler velocity is away from the observer and 
into the Sun. 
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Figure 13. A histogram of the duration of all SCB detections in all three events, 
minutes (dot-dashed line) and a median duration of 3.1 minutes (dashed line). 



The SCBs have a measured mean duration of 5.7 
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Figure 14. Derived statistics for all SCBs for the 6 May 2005 and 13 May 2005 events. All panels show the distance at which the SCB 
occurs versus time from the Ho flare peak. The color of each plotted point is representative of its relative intensity; the lighter points are 
higher intensity detections. Bottom: the same plots as the top two panels except with a Doppler filter applied - each point had to have a 
SCB Type I Doppler shift associated with the line center brightening. The dashed lines are fit to the two visually distinct groups in each 
event. 
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Figure 15. A plot showing the RHESSI non-thermal 25 - 50 keV intensity in black and the Ha SCB integrated intensity curve in gray 
for the 13 May 2005 flare. The vertical dashed line indicates the peak of the integrated SCB curve. 
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Figure 16. A cartoon model of the physical cause of SCBs. A dashed line marks the neutral line in this model. The central loops 
represent the emerging two ribbon flare arcade. The arches represent field lines connecting outside of the flare ribbons. The arrows on top 
of both the loops and arches show the direction trapped plasma flows as the flare begins to erupt and the stars suggest the location where 
the plasma is disrupted. Once perturbed, the trapped plasma streams down the loop lines and impacts the chromosphere causing the Ha 
intensity brightening. The different orientations of off-flare loops accounts for the different propagation speeds of SCBs observed. 
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Figure 17. A diagram of the physical dynamics occurring in a single SCB. Electrons and protons accelerated by magnetic reconnection 
further up in the corona come streaming down the flare loop. Since the mean free path of electrons is significantly less, the electrons deposit 
energy into the upper chromosphere. Unable to radiate the energy as heat, the chromosphere explosively responds by sending material 
back up the flux tube (chromospheric evaporation) and a reaction shock propagates towards the photosphere. This reaction shock explains 
the velocities observed in Type I SCBs. 



